Similarities between the Hubbard and Periodic Anderson Models 

at Finite Temperatures 
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The single band Hubbard and the two band Periodic Anderson Hamiltonians have traditionally 
been applied to rather different physical problems - the Mott transition and itinerant magnetism, 
and Kondo singlet formation and scattering off localized magnetic states, respectively. In this pa- 
per, we compare the magnetic and charge correlations, and spectral functions, of the two systems. 
We show quantitatively that they exhibit remarkably similar behavior, including a nearly identical 
topology of the finite temperature phase diagrams at half-filling. We address potential implications 
of this for theories of the rare earth "volume collapse" transition. 
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PACS numbers: 7f.27.+a 71.10.Fd 71.30.+h 

The single band Hubbard model (HM) and the two 
band Periodic Anderson model (PAM) have a long 
history in the study of many body effects in solids. 
While they have been proposed to describe complex d- 
and f-electron systems, respectively, the detailed atomic 
physics and orbital structure are nevertheless eliminated, 
leaving minimal Hamiltonians where the interplay be- 
tween electron kinetic energy, screened electron-electron 
repulsion, the Pauli principle, temperature, and electron 
density still give rise to a rich variety of phenomena. 

The key features of the half-filled HM are antiferro- 
magnetism and a Mott insulator-metal transition within 
the paramagnetic phase |l|,||. Originally ||, the Mott 
transition was associated with a change in the density 
of states (DOS) from widely separated upper and lower 
Hubbard bands to a broad one-peak structure, and an 
accompanying loss of local moments. The key features of 
the PAM are magnetism arising from an indirect 'RKKY' 
interaction of the local moments mediated by the con- 
duction electrons, and a DOS which, in the paramag- 
netic phase, has a three peak structure - upper and lower 
bands and an additional Abrikosov-Suhl resonance near 
the Fermi surface Q . Concomitant with the development 
of this resonance, local moments become screened by the 
conduction electrons, and an anomalous sharp feature in 
the free energy is observed S. During the past years, 
it has been established by dynamical mean-field theory 
(DMFT) @,§ and quantum Monte Carlo (QMC) simu- 
lations that the HM also has a three peak structure 
in its paramagnetic DOS, thereby blurring at least one 
of the features believed to distinguish the two models. 

It is the purpose of this paper to note even greater simi- 
larity between the two models at finite temperature, espe- 
cially if a more physical choice is taken for the interband 
hybridization in the PAM. Specifically, our new DMFT 
results show that the f-electrons of the PAM undergo a 
transition similar to the Mott transition of the HM, as is 



reflected by the f-electron spectral function, quasiparti- 
cle weight, and charge compressibility. Furthermore, the 
phase diagrams, at finite temperature, associated with 
the two models are topologically equivalent . 

These similarities are of particular interest in that the 
HM H and PAM (or its impurity approximation have 
been proposed as competing models for the "volume- 
collapse" transitions in the rare earth and actinidc met- 
als. These pressure-induced phase transitions are charac- 
terized by unusually large (9-15%) volume changes, and 
are believed to arise from f-electron correlation effects 
|lO|,[Ll]]. The present paper suggests that the predictions 
of the two models are more similar than has previously 
been thought, and that perceived incompatibilities may 
have arisen in part from differing and less rigorous ap- 
proximations used in the treatment of electron correla- 
tions, as for example in mean-field-like local-density func- 
tional theory |ll[] . 

The single band HM is given by, 
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Here e k = £g — 2i//[cos(fc x ) +cos(fcj,) + cos(fc z )] is the 
nearest-neighbor band dispersion, where the energy scale 
is set by t // = 1 , and Uf is the f-electron local interaction 
on every lattice site i. The two band PAM is, 

H = 4<A. + e{ flf ka + V k (d\j ka + h.c.)} 

ka 

+UfY,{nif1-\){nin-\), (2) 
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where Vfe denotes the hybridization between f and d elec- 
trons at momentum k. In the PAM, the d-electrons are 
chosen to have the same nearest-neighbor dispersion as 
the f-electrons of the HM (tdd = 1), while the f-electrons 
are chosen to be dispersionless, tff = 0. We study the 
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"symmetric" PAM, i.e., e 



d(f) 



and chemical potential 



|U = 0, which has the property that (n/) = (rid) = 1 
for all choices of Uf, tdd{ff), ^/dj an d temperature T. 
Finally, we choose a nearest-neighbor intersite hybridiza- 
tion Vk = —2tfd [cos(fc x ) + cos(fcy) + cos(k z )], but also 
have investigated the more canonical onsite hybridization 
Vk = VStfd, as well as an odd-parity nearest-neighbor 
form Vk = —itfd [sin(/c2;) + sin(fe y )+sin(fc z )]. The three 
choices have zero, constant, and maximal | Vk \ on the half- 
filled Fermi surface, respectively; and yield metallic, insu- 
lating, and semi-metallic /-projected densities of states, 
respectively, at U — and within a Fermi liquid phase. 
We favor intersite 14 due to the similarity of its metal- 
lic f-DOS to paramagnetic local density results for the 
rare earths and to the character of their f-electrons the 
high pressure phases. However, it is to be emphasized 
that the essential results of this paper are independent 
of these hybridization choices. 

Our calculational approaches are dynamical mean- 
field theory (DMFT) P,|q,[l2| and determinant quantum 
Monte Carlo (QMC) |13|jr4[ ], carried out for a three- 
dimensional simple-cubic lattice in the thermodynamic 
limit and for a 4 3 periodic lattice, respectively. DMFT 
can capture thermodynamic phase transitions while the 
quantitative agreement with QMC provides compelling 
evidence its accuracy. 

Within DMFT, both the f-electron part of the PAM 
and the HM map onto the same single site problem in 
which the Green function is calculated from a dynamical 
mean field S/(w) [§Ml5|. The difference between PAM 
and HM is the way e k and Vf. enter the self-consistency 
condition, i.e., the k integrated Dyson equation 
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where V k = and e{ — for the HM and PAM, re- 
spectively. In the DMFT equations the d-electrons of 
the PAM, which enter only quadratically, have been in- 
tegrated out yielding an effective f-electron problem with 
retarded potential 14 2 /(w — et). One expects the differ- 
ent potentials to be reflected in quite distinct physical 
properties of the two models. 

The spectral function A / (uj) of the f-electrons is shown 
in Fig.l. At small hybridizations tfd, the PAM spectrum 
consists of two peaks at ±{7/2 which corresponds to a 
phase with localized f-electrons. At larger tfd, an addi- 
tional Abrikosov-Suhl resonance develops at the Fermi 
energy which is not gapped for the PAM with inter- 
site hybridization. Thus, the f-electrons are itinerant. 
The striking similarity to the spectral function of the 
HM (inset) calls for a more quantitative comparison of 
the crossover between itinerant and localized f-electrons 
which is driven by an increase of Uf /tfd(ff)- Such a de- 
tailed comparison is the subject of this paper. 
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FIG. 1. The spectral function of the PAM at U f = 6 and 
T — 0.2 with an intersite hybridization of (a) tfd = 0.3 and 
(b) tfd = 1, respectively, as obtained by DMFT (solid line) 
and QMC (dashed line). With the exception of a small AF 
splitting of the central resonance found in finite-d simulations 
at tfd = 1 the agreement between DMFT and QMC is excel- 
lent. Insets: HM (DMFT) at U f = 6, T = 0.2, t ff = 0.3 and 
tff = 1, respectively. 

In the self energy T,f (Fig. 2a), the crossover shows 
up as an abrupt change of the behavior with tfd'- for 
large tfd Fermi liquid behavior with Iva£,f{ioj n ) oc uj n at 
low (Matsubara) frequencies occurs, while at small tfd 
an insulating-like behavior with diverging ImEj is ob- 
served. From the data one can calculate a quantita- 
tive measure for the crossover from itinerant to localized 
f-electrons: the quasiparticle weight of the Abrikosov- 
Suhl resonance Zf = (l — dTif(u})/du>\ LL ,=o)~ 1 (calculated 
via finite differential quotient (1 — ImEf(iu>g)/ uio) ) . 
Just as at the Mott transition of the HM 0,0|J^Q . 
Zf vanishes with decreasing tfd/Uf for the PAM. This 
behavior of the PAM is astonishingly similar to that seen 
in the HM and is insensitive to the choice of hybridization 
as is seen in Fig. 2b. 

While the important similarity between the two mod- 
els lies in the presence or absence of the central reso- 
nance as a whole (measured by Zf), the PAM hybridiza- 
tion can determine finer details such as whether this 
resonance may be split by a gap. Consider the charge 
compressibility k/ = dnf / 'd/j, shown in Fig. 2c. Nonzero 
Kf indicates a metallic phase, while Kf is exponentially 
small in 1/T within a gapped insulating phase. At small 
t 



fd(ff) 



all / spectral functions show insulating behavior 



as there are no central resonances lying between the two 



2 



Hubbard bands. These resonances form in all cases at 
larger however, they are split by insulating and 

semimetallic hybridization gaps for the onsite and odd- 
parity cases, respectively, as is evident in Fig. 2c. 
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FIG. 2. (a) (Matsubara) frequency dependence of the 
imaginary part of the self energy of the PAM with intersite 
hybridization. The behavior for the HM (not shown) is very 
similar, (b) Quasiparticle weight Zj as a function of t — tfd 
and t = tff for the PAM and HM, respectively (□ : PAM with 
intersite, A: PAM with onsite, and V: PAM with odd-parity 
intersite hybridization; ■: HM). (c) Electronic compressibil- 
ity Kf = drif/dfi as a function of t. All results are calculated 
by DMFT at Uf = 6 and T = 0.15 (A, V: T = 0.07). 



While the (unscreened) local moment 



((n T 



ny) ) shown in Fig. 3 decreases with increasing hybridiza- 
tion, it is still quite substantial through the HM and 
PAM crossovers of Fig. 2, marked by arrows in Fig. 3. This 
smooth behavior of the local moment from weak to strong 
coupling in both models is insufficiently appreciated in 
the discussions of the distinctions between the 'Mott' 
and 'Kondo' scenarios for the volume collapse transition, 
where mean field treatments of the correlations lead to 



artificially abrupt changes in the local moment. Despite 
the finite moments, screening effects lead to Pauli like 
low-T susceptibilities at large £/<z(//)- 




FIG. 3. The square of the local moment is shown as a 
function of t = tfd and t = tff for the PAM with intersite 
hybridization (open symbols) and HM (closed symbols), re- 
spectively, at Uf = 6 and T = 0.15 (circles: DMFT, squares: 
QMC) [T = 0.5 for the QMC HM to avoid antiferromagnetic 
ordering]. In neither case does the moment change remark- 
ably in the vicinity of the transitions of Fig. 2 (arrows). 

The T-Uf phase diagrams which result from consider- 
ation of the self-energy, quasiparticle weight, and mag- 
netic susceptibility are shown in Fig. 4. Both the HM 
and PAM have a low temperature phase of antiferromag- 
netic order with a maximum in the Neel temperature 
Tjv at intermediate Uf/W. At a slightly larger Uf/W, 
a second transition, i.e., the disappearance of the quasi- 
particle peak, is observed within the paramagnetic phase 
(if antiferromagnetism is frustrated). This is the Mott 
transition of the HM and, from the Kondo point of view, 
the crossing of the Kondo temperature Tx{tfd)- Apart 
from the temperature and hybridization scales, the phase 
diagrams of HM and PAM are remarkably similar at fi- 
nite T. It is an open question at present whether such 
similarities will persist down to zero temperature. 

The present paper has shown a remarkable similarity 
between the HM and the f-electron system of the PAM 
which is indeed - after integrating out the d-electrons - 
like a HM with a more complicated retarded potential. 
Novel numerical results on the spectral functions, charge 
compressibilities and local moments show that despite 
its more complicated potential, the PAM undergoes (at 
finite temperatures and half-filling) a crossover from itin- 
erant to localized f-electrons similar to the Mott transi- 
tion of the HM. While this similarity is closest for the 
PAM with intersite hybridization which describes metal- 
lic f-electrons at large tfd the vanishing of the central 
resonance does not depend on the hybridization choice. 
This crossover and the transition towards antiferromag- 
netic order yield phase diagrams which have the same 
topology for both models. Accordingly, the present pa- 
per suggests that the physics underlying the Mott f| 
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FIG. 4. T-Uf phase diagram of the Hubbard model with 
a semi elliptic DOS (filled symbols; reproduced from | pp| , ^8| ) 
and PAM with intersite hybridization (open symbols) within 
DMFT. Here W is the band width of the Bethe DOS and 
12 tfd for the PAM, respectively. Squares: Neel temperature; 
triangles: Mott transition/Kondo temperature (the dashed 
lines end at higher T as the crossover becomes too gradual in 
Uf/W to allow a precise location). 

(HM) and the Kondo volume collapse scenarios || (PAM 
or its impurity approximation) for the volume collapse 
transitions in f-electron metals may be similar and 

not incompatible. While the Mott transition is continu- 
ous at the temperatures considered, the rapid change in 
the correlation energy, when combined with the smooth 
contributions from other valence bands, might neverthe- 
less still result in a Maxwell construction giving the ther- 
modynamic instability manifested as the first-order vol- 
ume collapse in Ce pl| . 

In regard to local density calculations, orbitally- 
polarized and self-interaction-corrected applications p3] 
to the volume-collapse transitions resemble a mean field 
solution of the HM, where the transition coincides with a 
coalescence of the two Hubbard bands and a simultane- 
ous loss of local moment. The implication of the present 
paper is that regardless of the balance between f-f (HM) 
and f- valence hybridization (PAM), a correlated solution 
should show some Kondo-like attributes in the collapsed 
phase, including an analog of the three-peak local DOS, 
persistence of the local moment, and its screening as de- 
duced from the susceptibility. Local density theory does 
not show such Kondo-like features in the collapsed phase, 
however, it does nevertheless yield excellent structural- 
dependence of the total energy |23] . 

We acknowledge valuable discussions with N. Bliimer, P. 
van Dongen, M. Jarrell, and D. Vollhardt. Work was sup- 
ported in part by a DOE ASCI grant, the DA AD, the LLNL 
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